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BASIC NOTATIONS 


a = semi-major axis of the orbit of Jupiter 

a' = semi-major axis of the orbit of Saturn 



/j = secular part of / 

/p = periodic part of / 
i = orbital inclination of Jupiter 
i' = orbital inclination of Saturn 
= Gaussian constant 
m = mass of Jupiter 

m' = mass of Saturn 

n = mean motion of Jupiter 

n' = mean motion of Saturn 

^Pi Pi ~ Newcomb’s differential operator, associated with the expansion in powers of e 

^'pi Pi ~ Newcomb’s differential operator associated with the expansion in powers of e' 

r = radius vector of Jupiter 

r' = radius vector of Saturn 

1 f oos S 

R ; — = disturbing function for Jupiter 

P r'^ 

R' GOSS ' _ function for Saturn 

P 

1 ' 

Rq = — = direct part of the disturbing function for Jupiter and Saturn 


V 



^ ^ ~ indirect part of the disturbing function for Jupiter 

COS S 

/? C = — = indirect part of the disturbing function for Saturn 


S 

X 

x' 

X, 


= angle between r and r' 

= (X^,X2,X^) 

= (X[,X'2,X'^) 


= eiX 
= e'jV 

= e^K'^ 
= e^X'"i 


- Hansen coefficients 

P>Q 

7 ' -iyvyvy\) 

= mean longitude of Jupiter 

y*y = mean longitude of Saturn 

(In the later part of the exposition, these notations designate the canonical variables.) 
a 

a = — 


= -+a 
a 


= — a 
a 


(Beginning with Section 3, 7 designates 2 sin z/2.) 


= 2 sin 2 


7' = 2 sin^ 


VI 



r 


71 

y'l 

72 =^7e‘"'^ 

72 =^7'e’^'®' 

e = orbital eccentricity of Jupiter 

e' = orbital eccentricity of Saturn 



6 = longitude of the ascending node of the orbital plane of Jupiter 

d’ = longitude of the ascending node of the orbital plane of Saturn 

X = 

X' = 

M =k^(l+m) 

jLt' =k^(l+m') 

u = 2n - 5n' = mean motion of the critical argument 
7 T = longitude of the perigee of Jupiter 

7t' = longitude of the perigee of Saturn 

p = distance between Jupiter and Saturn 
X = = critical argument of the great inequality 


I 



w 

CO' 

_d_ 

dx' 

dy 

d 

dy' 


= argument of the perihelion of Jupiter 
= argument of the perihelion of Saturn 

= gradient operator relative to x 
= gradient operator relative to x' 

= gradient operator relative to y 
= gradient operator relative to y' 


viii 



THE INFLUENCE OF THE GREAT INEQUALITY 
ON THE SECULAR DISTURBING FUNCTION 
OF THE PLANETARY SYSTEM 


by 

Peter Musen 

Goddard Space Flight Center 


1. INTRODUCTION 

This paper discusses the influence of the great inequality, between Jupiter and Saturn, on the 
secular (long-period) disturbing function of the principal planets. Earlier results concerning this influ- 
ence were obtained by Hill more than 70 years ago. 

The great inequality is a near-resonance effect produced by commensurability close to 2:5 be- 
tween the mean motions of Jupiter and Saturn. 

The corresponding critical argument 2y^ - Sy'j in the trigonometrical expansion of the disturbing 
functiori has a period of approximately 900 yr. So long a period and close a commensurability pro- 
duce large perturbative effects in the elements, especially in the mean longitude of Jupiter, with ampli- 
tude nearly 1 200", and Saturn, with amplitude nearly 2900". Such periodic perturbations produce, in 
higher approximations, an appreciable effect on the secular disturbing function and on the secular 
behavior of the elements of the principal planets. 

Secular perturbations are the source of the long-period effects, with periods ranging from 5.7 X 10“^ 
to 2 X 10® yr. 

These perturbations are chiefly responsible for the behavior of the elements of the principal 
planets over millions of years. The amplitudes and mean motions of the arguments of these periodic 
terms determine the secular change in the longitudes of the perihelia and nodes. The amplitudes also 
provide information about the range of oscillation of the eccentricities and inclinations. Knowledge of 
the long-period perturbations of the motion of the principal planets is essential if we are to understand 
the secular behavior of the asteroidal ring, because these perturbations induce the forced long-period 
oscillations in asteroid elements. 

Since Lagrange, the classical way of treating secular changes in the motion of the principal planets 
has been based on the linearization of Poincard’s canonical elements or some of their approximations. 
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Only terms quadratic with respect to the eccentricities and inclinations are retained in the expansion 
of the secular disturbing functions. 

Modern computational technique permits terms of higher order relative to the eccentricities and 
inclinations to be included in the secular disturbing function and in the perturbations. Harzer (1895) 
in his classical work on secular perturbations solved the Unearized problem and made an attempt to 
include terms of higher order, but unfortunately this last part of his work was never completed. Only 
recently Anolik, Krasinsky, and Pius (1969) succeeded in developing the trigonometrical theory of 
direct secular perturbations (of rank zero) including terms of the fourth order with respect to the 
eccentricities and inclinations. 

Secular perturbations of rank 1 produced by the great inequality were omitted from their work. 

Hill (1897) tried to derive these perturbations from the Leverrier (1874) differential equations 
for the eccentricities and perihelia of Jupiter and Saturn. 

Brouwer and van Woerkom (1950) and Sharaf and Budnikova (1967) incorporated Hill’s contri- 
bution to the secular disturbing function in their linearized theories of secular perturbations of major 
planets. 

Brouwer and van Woerkom, however, assert the following “[Hill’s computations] cannot be con- 
sidered as definitive, because many of the coefficients were not derived in a rigorous manner. They 

were obtained as the means of the coefficients arising in four equations While in some cases the 

agreement was satisfactory, there were contradictions in other cases. These contradictions were re- 
solved by Hill in an empirical manner. ” 

In other words. Hill’s results need revision. We can fully understand the difficulties he encoun- 
tered, because at that time the theory of elimination of periodic terms from a Hamiltonian was not 
fully developed. 

Owing to the publications of von Zeipel (1916), Brouwer (1959), Hori (1966), and Deprit (1969) 
on the modernization of the method of Delaunay, we now possess an easy and exact algorithm for the 
elimination of short-period terms from a Hamiltonian. Even when the equations of motion do not 
have a canonical form, short-period terms can be eliminated directly from the differential equations by 
the method of Krylov and Bogolubov (Bogolubov and Mitropolsky, 1961). This method is described 
by Musen (1965) in a form convenient for application to astronomical problems. 

We may consider the great inequality as a short-period term compared with the secular effects. 

We can eliminate the great inequality either from the Hamiltonian or from the differential equations 
for the variation of elements. This elimination results in the introduction of secular effects of rank 1 
and class 1 /2 into the expression for the perturbations of the elements. 

HiH’s addition to the secular Hamiltonian contains terms of the fourth and sixth order with re- 
spect to the mean orbital eccentricities of Jupiter and Saturn. We can add those terms of the fourth 
order that depend on the inclinations (they are omitted in Hill’s exposition). The use of heliocentric 
elements of motion is the most convenient from an astronomical standpoint. 
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In our problem we can write the equations for variation of heliocentric elements in a canonical 
form for the whole system only if we retain the direct parts of the disturbing functions and neglect the 
indirect parts. This is because the indirect parts contain great-inequality terms of the fifth and higher 
orders and because they are different for Jupiter and Saturn. 

In order to retain the heliocentric elements and the canonicity we must neglect such terms. How- 
ever, we retain those terms of the third order that depend on the great-inequality argument, because 
they appear only in the direct parts of the disturbing function. Thus, in the secular disturbing func- 
tion we can retain terms of the second, fourth, sixth . . . order and of rank zero, also terms of the 
fourth order of rank 1 and class 1/2. These last terms are produced by the. elimination of the great- 
inequality terms of the third order mentioned above. 

If we decide to include in the disturbing functions terms of the sixth and higher orders in e, e', 7 , 
and 7 ' and to retain the canonicity of the equations of motion, then we shall use Jacobi’s reduction of 
the differential equations of planetary motion to the canonical form and use some modern version of 
Delaunay to eliminate the periodic terms. 

If the motion is referred to the Sun and we wish to retain terms of the sixth and higher orders, 
we can eliminate the periodic terms by the method of Krylov and Bogolubov. This choice is more 
convenient from an astronomical standpoint. 

Hill’s empirical adjustment of the coefficients stands between these two methods and evidently 
does not solve the problem. In our exposition we limit ourselves to terms of class 1 /2 and, therefore, 
neglect the secular effects produced by short-period terms in higher approximations. 

We found that Andoyer’s (1923) expansion of the planetary disturbing function suits our purpose 
very well, because he expands the function in powers of the inclinations relative to a fixed ecliptic. 

The transition from Andoyer’s expansion to Harzer’s expansion in the canonical elements is 
very easy. 

We found, in general, that many theories developed by Andoyer now sound quite modern; we 
believe that interest in his remarkable works should be revived. 

2. ANDOYER EXPANSION OF DISTURBING FUNCTION 

The Andoyer expansion of the direct part of the disturbing function is given by 


TWO 


where 


A{ =y xPixPAx'Ax'AN^^*> N'yV' 
S ^ 1 2 1 2 Pi,P2 P1.P2 
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(2a “ 1)*^ 


,q^\q^\q\\q'^ 


?1 + ?2 ■*■ ^1 ■‘‘ ^?2 = ^ > 




<^1 = (Ti - 7'i)72 " (72 " 72)7i - 7i72(7i " 7i)(72 " 7^ + ‘ ' ' , 


= (72 - 72>7i - (7i - 71)72 " 7271 (7i " 7l)(72 " 7l> + ' ' ‘ , 

«^1 = (7i - 7l)^ + 7i7l(7i - 7l)(72 “ 7^) + ' ' ' , 

and 

= (72 - 7l)^ + 727l(7i ~ 7l)(72 ~ 7l> + ' ' ’ • 

The definitions of Laplace coefficients and Newcomb operators given by Andoyer differ slightly 
from the classical ones. 

The term 

xfixf2xlPixlP2 

is of the order 


Pi + Pj ■*■ p'\ ^ ’ 

whereas 5^ is of the order 2/ in y and y'. 

The Laplace coefficients and their derivatives in the Andoyer expansion are successively com- 
puted from the following equations; 


-IT 1 2 

I ( 1 - cos^ cos^” ij/ d\J/ 


^n+1 > 




^P+1 = —DbP , 

n py n 


Z)26P = (K±p)26P-h4p26P+l , 


and 


DbP*^ =--{-DtP„ +D'^b^ , 

P7\7 / 

= (« ± p)'^D^bP^ + Ap^D^bP^li . 


The differential operators and are computed from the following recursive equations 


(Andoyer, 1923): 


and 


where 


We have also 


and 


.P2 


,P2 


m ^ V 

- 3(a + p, - pp 2 ^ 

Po^n n = r20+4-^Vr‘„ -1 

+ (-a+^-Z))Af^-2 
\ 2 / Pl>P2”^ 

4 3(.4p, -p,) I . 


m = min (ppPj) • 


Pl»P2 P2>Pl 


wjo = > . 


Wi" . (+Z)) = W" „(-B) 

PI.P2 PI.P2 


(1) 


( 2 ) 
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C. Hipkins of the Analytical Mechanics Associates prepared a program, using an analytical lan- 
guage, for the computation of operators N and N' in rational-fractions arithmetic. In his program, N 
and N’ can be obtained up to any desired order in the form of polynomials in a and D. The use of 
both equations, (1) and (2), provides a sharp check of the whole calculation. Hipkins found a misprint 
in (2) as it is given by Andoyer. Instead of 


prp>P2~p 

under the summation sign in the book, one should read 

prP’P2~P 

Such a program paves the way for the production of analytical planetary theory on electronic machines. 

3. TERMS ASSOCIATED WITH THE GREAT INEQUALITY 

If we want to extract from the disturbing function those terms having arguments that are multi- 
ples of the given argument iy-^ + or, what is the same thing, to extract those terms having the 
factors 


X-"' , X'^'' , X~' , and X'"*' , 

from Andoyer ’s expansion up to the order P in the eccentricities and inclinations, then Pj, pj* Pp 
p '2 must be selected in such a manner as to satisfy the following conditions: 

Pi +P2 +Pi +P2 


and 

Evidently the inequality 
or the equivalent inequality 


Pj -P2 + s+j = ik , 
p[-P2-s+j = i'k, 
5,/,k = 0,±l,±2, ... 


IPi -Pi\ + Ip'i -P2I 


\ik - 5 - ;| + \i'k + s~ j\^P 


must also be satisfied. The use of the last inequality facilitates the search for the prospective candi- 
dates for k, s, /, Pj, P2, Pj, and p^. In our case, 

i = 2 , i' = - 5 , 
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and Table 1 shows the admissible values of the indices and exponents. To the required accuracy we 
obtain 


+ V-4 VSa) + V-5 

+ iAllBt\ + A:\BII) + (AllBl\ ^Al\B-\) , 

or, if we substitute 

Xj = , Xj = X"^e2 , Xj = X'^^^e'j , Xj = X'^^e^ 

and take, with the necessary accuracy, 

B^,^ = V2o\bl!\ 

and 

s;‘ = . 

we obtain for the direct part of the disturbing function associated with the great inequality the follow- 
ing expression; 

= (X'^X'-^^e^ + 

+ (X^^X'+^e^e'j + X+^X'-^e^e^Pj 


Table 1 —Admissible values of exponents and subscripts. 


/ 

s 

k 

Pi 

■ 

Pi 

P'l 

p'l 

0 

-2 

-1 

0 

0 

3 

0 

0 

+2 

+1 

0 

0 

0 

3 

0 

-3 

-1 

1 

0 

2 

0 

0 

+3 

+1 

0 

1 

0 

2 

0 

-4 

-1 

2 

0 

1 

0 

0 

+4 

+1 

0 

2 

0 

1 

0 

-5 

-1 

3 

0 

0 

0 

0 

+5 

+1 

0 

3 

0 

0 

+1 

-3 

-1 

0 

0 

1 

0 

-1 

+3 

+1 

0 

0 

0 

1 

+1 

-4 

-1 

1 

0 

0 

0 

-1 

+4 

+1 

0 

1 

0 

0 
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+ (X"2X'-"^e;3 + X-"2X'"5e^3)p^ 

+ (X"2x'+5g^(j'^ + X‘*’2X'"^e2ap0Q 
+ (X"2X'-"^e;a'i + , 

(3) 

where 




13 2_i^3\^l/2 

\ 48 24 4 ^ 6 /^5 ’ 

(4) 


^''2,0^^ 1,0 ^4 




(5) 


‘"2^N-,lN'^bf 



/_2585_587 43 j_I 

V 16 8^4^ 2^/^3 ’ 

(6) 





= /'+2«5 6U 15 3\ „2 

V 48 ^ 24 ^ 4 ^ 6/2 > 

(7) 


/O = +-17V~3a^'+5a3/2 
^0 2 '3,0 ^^4 



= (-^-^>f. 

(8) 

and 

O = +l-A^+2;u^f+4L3/2 
^1 2 0,0'''^ 1,0 ^3 



= ( 4 ^ 40^6,3/7 , 

(9) 

The first four terms in 

(3) and the coefficients in (4)— (7) are given by Andoyer. 
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r 


For the indirect parts of the disturbing functions we have (Andoyer, 1 923) 




2 PvP 2 Pl>P 2 J 


and similarly 


+ a;x"^x'‘^x-2--izVgil . 

2 Pl>P2 Pl’P2 J 

The great-inequality term of the lowest order in i? j is 

+ (1 + a2)X-lX'^lx2X;^Z+j--lZ4-2. + l] , 


and in the lowest term is 

All these terms produce terms of the sixth order in the secular disturbing function; and they can 
be omitted if we decide to retain the canonicity of the differential equations for the heliocentric 

elements. 

r 

There will be a large number of sixth-order terms. This usually means that they have a small 
1 weight in the combined secular effects. 

4. ELIMINATION OF THE GREAT INEQUALITY 
FROM THE HAMILTONIAN 

In performing the elimination of the great inequality we can consider only Jupiter and Saturn. 
For the Hamiltonian we can take 
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m^l m'lx' fmm' 

and neglect terms above the third order in the expansion of F. 

The canonical elements of Delaunay in our case are 

L = my/Jm = mna^ / = mean anomaly of Jupiter 

G = L^\- g=oj 

H = G cos / h = 6 

L' = m’y/ix’a' = m'n'a'^ V = mean anomaly of Saturn 


G' = L'^\ -e'2 g' = co' 


H' = G' cos /' 


h' = 6' 


Dividing F into purely secular and purely periodic parts, we have 


where 




3 2 t?, f'l 

p = 1 

° 2I2 2I'3 


Pxs ^ f 1^2 + ^'l^2 - (^1^2 + 

+ {e\e\ + o\e’2 + 

+ + cr'^Cj + 4^0^^ + Db]^^ 

- + e2e\)(^b^f +-|*p -Dbpj 

+ 916^626162 ■'■°'2)(^1^2 ■*■ ^0^^ 


( 10 ) 


10 


and 


where 


and 


■f + <^2«1^2 + ‘^1^2 ■*■ 

+ (aje^e^ + c^2^2^i)(*P -f ^p)|. 


= (7-X+2x'- 5 + oX-2\'+5) ^ 

^ / /Tr/Tf ^ 


( 11 ) 

( 12 ) 


= ^ 0^2 + ^ 1 ^ 2^2 + ^ 2 ^ 2 ^ 2 " + + ( Qo ^2 + QAlK - 

a = Poe3 +p,efe[ + P2e,e[2 + P^e[^ + (Q^e, + Q^e[)a[ . 

At the beginning we set up elimination of the great inequality, using the following canonical 
elements: 


(12') 


x^= L 
X2= L - G 
X2 = G-H 

jKj = 1+ g + h 

y2 = -S~h=--n 
y2=-h = -e 


x[ = V 
x\=L'-G' 
x\ = G'-H' 
y'l = I' +g' +h' 

y' 2 =-g'-h' = -Tr' 

y',=-h=-e' 


At a later stage we switch to the canonical elements of Harzer and Poincard, because of their intimate 
connections with Andoyer’s expansion. Harzer’s elements are 

x\=L' = g'^ 


yi=yi 

\b 

= 2gsinye 
^2 = 2g sin y 


y'i=y'i 

\b^ . , 

?'i = 2g' sin -y e~"’ 

\b' • / 

?2 -2g' sin-y 
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where we set 


Pj = 2gy/coT^ sin 



Pj = 2g'\/cos^ sin 



g = a-i/mn g' = a's/m'n' 

e = sin \]j e' - sin ip' 

In eliminating the great inequality we do not go beyond terms quadratic in fmm’ and thus must 
use very simple differential operators associated with the process of elimination. (We are not particu- 
larly interested in the inversion process of expressing the old elements in terms of the final ones.) 

For all these reasons we found it convenient to use von Zeipel’s method in its classical form to 
eliminate the great-inequality terms. In addition, when using this method we need compute at each 
step only one half of each Poissonian bracket to obtain the unknown secular terms. When using other 
methods we must compute full Poissonian brackets. The author is indebted to Dr. Garfinkel* for 
pointing out both these facts. 


and 


We make the canonical transformation 

d5(x*,y;x'*,y') 








, X. X. T 

, y'*=y'i + 


dS{x*,y,x’*,y') 

dy; 

dS(x*,y;x'*,y') 

dx;* 


(13) 


where 


/= 1,2,3, 


in such a manner that the new Hamiltonian F* does not contain the angles and yj*. 

The essence of von Zeipel’s method can be described in a contracted form as follows: changing 
notations in (13) from x*, x'* to x, x', we can write the relation 

F(x,y;x\y') = F*(x*,y*;x'*,y'*) 


in the form (von Zeipel, 1916) 

dSix,y;x',y') 


FU + 




d5'(x,y;x',y')\ 
a? ) 


= F*(x,i 


dS(x,y;x',y') 35'(x,y ;x',y')\ 

^ a? ) 


(14) 


*Garfinkel, B,: 1969, Private Communication. 
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I 


Introducing the Taylor operators 


and 


T= exp 


\9y 9x dy' dx') 

^ ^^^\dxdy dx'dyj’ 


we can rewrite (14) as 


TF(x,y;x\y)=T^F*(x,y;x\y^) . 


(15) 


Assuming the expansions 


5 = jSj + *S'2 + * * * , 


and 


F = Fq+F^+F 2 + --- , 
F* =F*+F\+Fl + --- , 


in powers of fmm\ we can obtain the expansions 


r = Tq + Tj + 7^2 + * • • 

and 

r* = r* + r* + r* + • • • , 


in terms of Faa-de-Bruno operators (1855) Tj and 7?. The operators T. and 7? are polynomials in the 
operators 


and 


i dy dx dy' dx' 


(16) 


^ dx dy dx^ dy' ’ 


(17) 


respectively. From (15) we derive the system of partial differential equations 


n n 

(n = 0,1,2,...) 

*=0 fc=0 

for the determination of 5^, S 2 , .... The values of F^, F^, F^ ■■■ must be determined in such a manner 
that S 2 , ... do not contain any secular terms. Further details on this approach to von Zeipel’s 
method and the general equations for operators T., Tf (/ = 0, 1, 2, ...) are given by Musen (1965). 
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In particular, we have 


and 


8,F^+F^=8lF*+Fl, 

(«2 + + ^2 = («l + + W + P': 


( 18 ) 


Substituting (10) for Fq and substituting + F^^ for Fj [where F^^ and Pip are given by ( 1 1 ) and 
(12), respectively] and taking (16) and (17) into consideration, we obtain 


and 


as, 

-^-^^Pis^Pip=Pl 


as, 1 / a^Fo a^F, 

4t^+25 

9X 2\ dxj 

as 


^%\(^sy/ds^ a ^as^ 3 

a^ivV^V Wi 9^2 9^2 


9y 


i_i_ A A+^ a Y 

3 axj ay'^ a^j ay^ 9x2 ayj dxy '■ 


lp+^U> 


'95, a 


as. 


^ ^ + J- +“■ -8-W. +F. 

a^2 ay2 axj ayj ax^ ay^ axj ay^y ^ 2 • 


as, 


(19) 


From (18) we have 


F*=F 


Is 


(20) 


and, making use of (12), we have 


(^X^2x-5 + ,jX-2 V+5) 

9x 


Consequently, after the integration, we have 

fmm' 

S, =- 

ivy/aa' 


(rX^2x'-5 _ aX'2X'+5) . 


( 21 ) 


Because F^^ and F^ are purely secular and Sj is purely periodic, the only contribution to in (19) 
comes from the terms 


14 



and 


Thus, 


^ hS, ^ dS, dF,^ ^ dS, dF,^ 
3>^i 9xj ^ by 2 9^2 ^ dy^ dx^ 9;^'^ 9x'^ dy'2 9xj 9;;'3 9x^ 


2\f 9x2 25^^,2j\^3^j • 


1 I 92Fo 92f, 


92F()\/95 iY| 

^A^/ i 


,9>»3 9xj 9y2 9x2 


The term 


!ii?f!£\,!£i!£k+?^i?£i£'U 

^ 9>>3 9X3 ^ 9yj 9x'j 9^2 9x2 ^y'3 ^^'3 /is 

is of the sixth order in e, e', 7, 7\ and thus can be neglected. Substituting (12) and (21) in the last 
equation, we obtain 






') 


9(t, a) 9(r, a) 9(t, a) 9(r, a) 


|9(;^2’-^2) 3 (^ 3 , X 3 ) 3 (^ 2 ’ x^) 3 ( 32 ^, x' 3 ) 

Transforming the last equation to Harzer canonical variables, we obtain 


F* = —I 
^2 


2 / fmm'\ 


9(r, a) 9(t, a) 9(r, o) 9(r, a) 

+ + ^ —7 +- 


9(^1, ^2) 9(PpP2> 9 (?'i>^2) 


( 22 ) 


With sufficient accuracy we can set 


= 2g€3 

ri = 2g'e; 

^2 = 2^62 

r2 = 2g'e' 

(N 

11 

s£T 

P[ = 2g'7l 

P2 = 2^72 

P'2 = 2g’y’2 


Substituting these values in (22) and taking 
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9r _ da_ _ _ dr^ _ do^ _ Zr_ _ da 

de^ de2 dy^ dj2 9^2 ^'^2 


into consideration, we obtain 


F*=- 



da ^ 8g \ 

de^ 372 


We set 

and 

Taking 

and 

into consideration, we obtain 
and 


+ 1 / ^ ^ ^ V 

^m'«'a'2^3e^ 3e'i dy’2dy[j ' 


fmm' 1 
2vaa' mna^ 

fmm' 1 
2vm' m'n'a''^ ' 


/(I + m) = rP'a^ 
/(I +w') = «'V3 


m'n 

2v{ \ +m)a’ 
mn' 

2v(l +m')a ' 


Substituting the values of r and a given by (12') into (23), we obtain 


= fmm'^^9PlK + PlK')ejel + (^’ 2 ^ 0 '^ 

+ (PjK + 9PlK')€'^^e'2^ + (SPqQqK + 2i/’i/c')(e^a' + ej-a^ 

+ 3P^(3PqK + P2K'){e\e2e'2 + ^ 1 ^ 2 ^!) 

+ 2P2^P\*<- ■*■ 3i’3K')(e2e'j^e2 H-e^e'^e^^) 

+ 4(/’jK +P^K')eje2eje2 + 3(PqP 2K + P^P2ic')(e^e'2^ + 

+ 2(P^QqK +i’2Qi'c')(e2^2°''i + 


(23) 
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- 4(K + K')6o^ie2 ' ^2^ 


- 4(K + k')Go2i(^1^2 ■*■ ■*■ ^2^ 

- 4(/c + K')Qje^e '2 • %(cTj + 02 )^ ■ (24) 


Taking into account (11), (20), and (24), and omitting useless terms, we obtain the transformed secu- 
lar Hamiltonian associated with Jupiter and Saturn: 

f* = fi+f*2 


fmm'^^ie^e2 + e'jC^ + + ViUj) +5^(ej62 + e2e'j) 

+ B^e\el +B^{e’2^a\ + e[^o’2) + B^e[^e’2^ + <^2^l> 

+ B2{e\€2€'2 + e^ele\) + B^{e^e\e’^ + e2e'j2e^) + B^e^e2e\e'2 


where 


+ 5joeie2 • V2(aj + 02 ) +5n^'i^2 ' '/^("i + ^2> +^i2‘"i^2 
+ 5j4(e2e^2 + + + eie'202 + ^2^201 +eie;a'2) 

+ 5 ig(aieje' + 02626^) + • >/2(aj + 02) + B,^{a\ + a^)] , 
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+1 > 

Bg = +-^bf + (4P2« + 4/’2k0 , 

B^Q = +-^b^J^-A{K + K')Ql, 
a y/a 

B^^=+-^b^J^-A{K + K')Q\, 

U y/ CL 

5j7 = +-^ bl>^ - 4(k + k’)QI , 
ay/a 

and 

^18 = +l^*f • 

Sa'^/a 2 

All elements in F* and F^ are in fact the “asterisk elements.” We change the notations and omit the 
asterisks. 


With the elements given in the work of Anolik et al. (1969), we obtain (using C. Hipkins’s pro- 


gram) 


?o =-3.42811 , 
Pj = + 17.1336 , 


P, =-28.3531 , 


P, = + 15.4786 , 


I I I I I 


I III II II mil 


I mil 


00 = - 1.95863 , 

01 =+3.76733 , 

F* 

jr^ = +0.924$99e\el 

+ 0.542825(6^2(1/ +e'^) 

+ 6.823706^26^2 

+ 0A9990l(a[el + a'^ej) 

- 4.59557(616262 + e ^ ele [) 

- 8.32l65(e^e[e'^ + e2ej26p 
+ 10.148561626162 

+ 2.51217(626^2 +g2g/2) 

- 0.662 104 ( 62620 '! + ei6'ia2) 

+ 0.0432945(j'ia^ 

- 0.0592693 • + a^)e^e2 

+ 0.1 14001 • ^ 2(01 + CT2)(eie2 ^ 2 ^ 1 ^ 
-0.219276 • 1 / 2(01 + 02)6'i6^ • 


5. CONCLUSION 

In order to meet the requirements of modern Celestial Mechanics in improving the theory of sec- 
ular perturbations of the principal planets, it is necessary to include the influence of the great inequality. 

The contribution of the great inequality to the secular disturbing function is considerable. A 
comparison between and Fj shows that the great inequality contributes the largest portion of the 
coefficients of many fourth-order terms in e, e', y, and 7 '. The part F 2 produced by the great inequal- 
ity can be added to the disturbing function given by Anolik et al. (1969). Better still, we can preserve 
the homogeneity of the theory and recompute the secular perturbations of principal planets from the 
very start, using the Andoyer’s symbolism. 
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We have already pointed out that the canonicity of the differential equations and the heliocentric 
elements can be jointly preserved only if we do not go beyond the fourth order in e, e', 7 , and 7 ' in the 
secular disturbing function. The use of heliocentric elements is important from an astronomical stand- 
point; if we wish to retain them and at the same time include the terms of sixth and possibly higher 
orders, we can use the method of Krylov and Bogolubov for eliminating periodic terms. 

We may expect, in general, that the Krylov-Bogolubov method will occupy an important place in 
developing the analytical theories of heliocentric planetary motions. 
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